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This paper focuses on a photovoltaic system for pumping water. The control 
strategy for this water pumping system is based on Takagi-Sugeno type fuzzy 
supervisors and sliding mode controller. The first generates the maximum 
power point current under varying climatic condition whereas the second 
allows tracking the reference signal produced by the fuzzy supervisor. The 
system includes a photovoltaic generator (PVG) followed by a DC-DC 
Converter, DC bus, an AC/DC inverter which is connected to the induction 
motor. This latter is coupled with a centrifuge pump. The induction motor is 
driven based on field-oriented control strategy. The Takagi-Sugeno type fuzzy 
supervisor predicts, depending on the variations of climatic variables such as 
irradiation and temperature, the optimum operating point for the photovoltaic 
source. The simulation results show the effectiveness of the proposed 
approach in transient and stationary regimes for different values of climatic 
variables. 
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1. INTRODUCTION 

In these last decades, the use of photovoltaic energy has grown rapidly due to its direct obtainability 
without causing pollution or noise. Referring to the studies that have been made by experts, solar energy is 
abundant in southern Tunisia, the average sunshine rate per day is around 7 kWh/m? in addition to 
3200 hours of sunshine per year [1]. All of these have prompted leaders to think about “photovoltaic water”. 
In Tunisia, the main reason delaying the use of photovoltaic technology for pumping water is a funding 
problem. However, an interesting program has been developed to exploit natural resources in particular the use 
of a photovoltaic generator for pumping water. The surface water resources are very limited, but deep-water 
resources are much more important. In addition, the majority of the population is settled around small villages 
with limited access to conventional energy sources; their water needs are relatively low and correspond to a 
moderate water consumption satisfying a certain household use. In reality, most type programs are supported 
by technical cooperation projects, in particular with Germany. For example, a photovoltaic micro-plant for 
pumping water for the irrigation of a palm grove of around 105 hectares has been installed in the Kebili region 
[2]. The project aims to install 14 photovoltaic pumps sponsored by the German Cooperation Agency (GIZ). 
In the literature, several works have been developed concerning the electric drive system best suited to water 
pumping applications with photovoltaic technology in the Saharan regions [3], [4]. 
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Various works have been interested to PV water pumping system [5]-[7]. In this type of process, 
different kinds of motors have been used such as DC motors [8], brushless DC motor [9] and AC motors [10], 
[11], the choice depends on some factors such as the efficiency, the reliability and the price. However, water- 
pumping systems based on AC motors, particularly induction motors (IM), are more attractive because this 
type of motor is more robust than other types. They have high efficiency, easy maintenance, and low cost [3], 
[4]. Field oriented control (FOC) and direct torque control (DTC) are the most used control strategies for 
induction motor drives [12]. It is well known that FOC makes it possible to control the IM in a similar way to 
that of the direct current machine. 

The performances of a PV micro-plant system depend on irradiation, temperature and type of load. 
However, maximum power point tracking (MPPT) has been used to maximize the extraction of power from 
the PV panels. In this context, several works are interested to the MPPT algorithms. There exist three big 
categories of MPPT algorithms, which are known as conventional methods, artificial intelligence and hybrid 
methods. Indeed, Mohammed et al. [13] presented perturb and observe (P&O) algorithm, the authors show that 
this method is simple to implement but the operating point remains oscillating around MPP and the operating 
point is not able to follow the maximum power point when there is a fast variation in environmental conditions. 
Kish et al. [14] presented the incremental conductance algorithm, this latest is simple for implementation and 
it converge rapidly. It is more efficient than P&O algorithm but the operating point oscillate around MPP. 
Chalok et al. [15] proposed an algorithm based on fuzzy logic to reduce the oscillations of the operating point 
but the implementation of the algorithm remains relatively complex. Idrissi et al. [16] neural network has been 
applied to track the maximum power point this algorithm provides rapid tracking and reduces oscillation but 
its complex for implementation. 

The aim of this paper is to develop a new MPPT algorithm for an efficient PV water pumping system. 
The proposed algorithm is based on the combination of Takagi-Sugeno type fuzzy supervisor with sliding mode 
controller. From the climatic variables, essentially irradiation and temperature, the fuzzy supervisor is able to 
predict the coordinate of optimum operating point of the photovoltaic source as well as the maximum converted 
power. This power is used as a reference for a sliding mode controller to extract the maximum power from the 
panel and generate the duty cycle to the boost converter. The proposed controller guarantees good performance 
less oscillation and fast convergence to MPP. The power computed by the fuzzy supervisor is exploited to predict 
the reference speed of the induction motor, in order to control the flow rate of the pump. 

The paper is organized as follows: Section 2 presents the block diagram description of the photovoltaic 
water pumping system; In the first part of section 3, the modeling of the photovoltaic panel is given; An MPPT 
algorithm, which results from the combination of sliding mode control with fuzzy systems, is given in the 
second part of section 3; In section 4, the different elements of the PV water pumping have been presented 
such as the voltage inverter, the induction machine operating under IFOC and the centrifugal pump; and In 
section 5, simulation results are given to illustrate the performance of the T-S fuzzy supervisor and the proposed 
approach. The paper is enclosed by concluding remarks. 


2. DESCRIPTION OF THE PUMPING SYSTEM 

The photovoltaic pumping station essentially comprises a source of photovoltaic energy, a DC/DC 
converter, followed by a DC bus connected to an inverter. The latter feeds an induction machine coupled to a 
water pump. An MPPT algorithm based on fuzzy logic and sliding mode is used to drive the DC/DC converter 
in order to extract the maximum power from the photovoltaic source. The induction machine is controlled 
according to the field-oriented control theory so that it behaves like a DC machine. The block diagram of the 
pumping process is shown by the Figure 1. 
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Figure 1. General diagram for a PV pumping system 
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3. PHOTOVOLTAIC PANEL MODEL AND MPPT CONTROLLER 

This section includes two sub-sections. The first concerns the modelling of photovoltaic panel and 
DC/DC boost converter whereas, the second is reserved to the description of MPPT controller and fuzzy 
supervisor. The MPPT controller is basis on sliding mode and fuzzy system. The stability of system is verified 
by Lyapunov approach. Then the fuzzy supervisor is described to predict the optimal values of the main 
quantities of photovoltaic panel. 


3.1. Modelling of photovoltaic panel and DC/DC converter 

It is well known that the photovoltaic panel is obtained by a series-parallel assembly of a set of cells. 
Several mathematical models describe the evolution of the photovoltaic cell as a function of climatic 
parameters, mainly temperature and irradiation [17]. In this work, we adopt the model given by Figure 2. 


Figure 2. Equivalent electric model for a photovoltaic panel 


The model of the photovoltaic panel is governed by (1)-(3). 


V+RsI 
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The current produced by the photovoltaic panel changes according to the temperature and the irradiation in 
accordance with the expression described by (4): 


Inn = (Ipnn + KAT) a (4) 


Tynan represents the nominal current supplied by the PV panel under standard climatic conditions of irradiation 
and temperature (G = 1000 w/m? and T = 25 °C). 
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Where: Jp is a reverse saturation current. 


Vice “_ log (A) (6) 
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V. = n,“Tlog (“= ") (7) 
Where, V,. represents the open circuit voltage and J,. is the short circuit current. Figures 3 and 4 respectively 
display the evolution of the power characteristics as a function of the panel voltage following variations in 
irradiation and in temperature. 

Before approaching the MPPT algorithm, it is crucial to introduce the modelling of the DC/DC 
converter. In fact, the main function of the DC/DC converter is to set the operating point by acting on the duty 
cycle “Z , the variation range of which is limited between 0 and 1 (0 < uw < 0C,,). In reality, there are several 


types of DC/DC converter such as boost, buck, buck-boost. The choice of converter is made according to the 
application. In this paper, we use a DC/DC Boost converter of which the structure is given by Figure 5. The 
mathematical model of the DC/DC boost converter is described by (8), 
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where 7,,L and Czrepresent respectively the resistance of the inductance, the inductance and the output 
capacitor. 
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Figure 5. Structure of DC/DC boost converter 


3.2. MPPT controller based on sliding mode and fuzzy supervisor 
To extract the maximum power from PV panel, we propose a new MPPT algorithm that combines 
sliding mode and fuzzy logic techniques as displayed in Figure 6. The T-S Fuzzy supervisor provides optimum 


panel voltage and current, the sliding mode approach is employed to truck the reference current and provides 
the duty cycle, which will be transformed into a PWM signal. 
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Figure 6. Block diagram of the MPPT algorithm 


3.2.1. Sliding mode control 


To extract the maximum power from the PV panel, we use a controller based on sliding mode to force 
the current to follow its reference value, which is computed by fuzzy supervisor. The output of controller 
generates the duty cycle uw to the boost converter. However, the sliding surface is defined by (9). 


S=e,+Gf edt (9) 
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Where, G is a positive constant, e; = Iypp — I,, Iyppis the reference current at the output of the fuzzy logic 
supervisor. The derivative of the sliding surface is given by (10) and (11). 


$=Ge,+ é, (10) 
$= Ge, +*1,(Yyy- (1-H) (11) 


Basis in the mode theories, the control law includes two terms. The first term is known as equivalent control 
law which is computed froms = 0. Whereas the second term is known as switching control law which is 
expressed as: Ksign (S). 


H = Heg t+ Hs (12) 


The expression of the equivalent control law, eg is computed from s = 0. 


§=Ge,+ 11. Vpv —(1-p2,)%) =0 (13) 
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To check the stability of the system, we use Lyapunov's theory. The candidate Lyapunov function is chosen as: 
V =-s? (16) 


however, V = SS 
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basis on (14) the last equation becomes as: 

V = S[-Ksign (S)] (18) 

V =—-K|S| <0 (19) 
therefore, the stability of the system is guaranteed. 
3.2.2. Fuzzy logic supervisor 

The fuzzy systems are considered as universal approximators. They can consistently estimate all 
continuous functions defined in compact domains. Fuzzy techniques are able to convert nonlinear function into 
amalgamation of linear functions. The fuzzy system is described by an assembly of fuzzy rules whose structure 
is as follows [18], [19]: 

if Z,is M,,and Z,is Mj, and Z,is Min then y, = yi,y2 = y) and y3 = ys (20) 


where, {Mi ;} describes the fuzzy sets, Z,(t) .....Z,(t) are the premise variables. y(t) € R»; c is the number 
of fuzzy rules. Each premise variable derives from a specific bounded domain such as: Z;(t)€[Z; minZi max]: 

For each ruleR;, a weight w; (z(t))is given which depends on score of membership function of 
premise variablesz,(t) in fuzzy sets Mj;. 


wi(z@) =| [my (©); wi(2@) > 0; fori =1,..6; 
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Mi; (z)(¢)) is the score of membership of z;(€)to the fuzzy setM;;. 


h,(z(t)) = ee 0 his) <tr Sse 


The output of the fuzzy system is expressed as (21). 
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In this context, we have exploited a T-S type fuzzy supervisor given in Figure 7. The latter allows to predict in 
function of climatic variables, such as irradiation and temperature, the optimal values of the main quantities 
such as voltageV,,y_,ef, Current! ypp and power p,e¢. Figure 8 shows the membership function of irradiation as 
displayed which has two subsets M1; and Mj2. The membership function for temperature as shown in Figure 9 
is chosen with tow subsets M2,;and M22. 
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Figure 7. T-S fuzzy supervisor 
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However, the T-S Fuzzy model of PV panel is well defined by the four following fuzzy rules: 
Vii = Vuppi¥ai = ImppiY3i = Puri (for i = 1 to 4) 
Plant rules: 
Rule 1: if Z,is My;,and Zzis M2, then Vypp = Vupp1Jupp = Imppiand P = Pypi 
Rule 2: if Z,is M;,and Zzis M2, then Vypp = Vypp2,tupp = Impp2and P = Pyp2 
Rule 3: if Z,is M;,and Zzis M2, then Vypp = Vupp3,Iupp = Impp2and P = Pyp2 
Rule 3: if Z,is My,and Zzis M22 then Vypp = Vuppalupp = Imppgand P = Pypa 
Were, 
Mya (22(t))= pin 


Zimax(t)-Z1(0) 
Zimax—Zimin 


22(0)-Z2min 
Z2max-Z2min 


_Z2max()-Z2() 
M22(z2 CO) sare an 


M21 (z2(t))= 


My2(z1(t))= 


So, the optimum voltage, current and power for a photovoltaic panel can be estimated by the following: 


Vurp = diar hi (20) Vapi (22) 
Iupp = Dias hi (z(t) )Iuppi (23) 
Prax = Pup = Dinrhi (z(t)) Puri (24) 
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4. EQUIPMENT OF THE PUMP THROUGH THE DC BUS 

After the DC/DC boost converter, the pumping chain comprises a DC bus, which feeds an induction 
motor through a three-phase inverter. The objective is to control the flow rate of the pump, which is a function 
of the rotational speed of the induction motor. This leads us to control the rotational speed of the motor. In this 
context, the field-oriented control technique is used to overcome the difficulties presented by the nonlinearities 
of the induction motor. 


4.1. Voltage source inverter 


The structure of the inverter, which consists of six IGBT switches, is shown in Figure 10. It is 
controlled by analogue PWM-type functions. The mathematical model of the inverter is given by (25). 
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= =1. 2. 2185 
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Where S;, 53,53, 51,52 and S3 are the IGBT control signals and Ug: is the input DC voltage. 
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Figure 10. Structure of the inverter 


4.2. Induction Motor model and field-oriented control 
The state model of the induction machine in the (d, q) rotating frame is given by the following [20]: 
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The expression of the electromagnetic torque is given by (27). 
_ 3PMsr 
Com = 5 (X%2.X3 — X4X4) (27) 


2 Lr 


The reference state model of the induction machine is based on the field-oriented control. It is well known that 
the principle of field control for the induction machine consists of keeping the rotor flux constant along the 
direct axis and cancelling the component of the flux along the quadrature axis [21], [22]. 


Oa, = cte (28) 


Dqr =0 (29) 
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The reference state model of the induction machine under the effect of field-oriented control becomes 
as (30). 


: k 1 
Xp = YX, twox.+—%x%3+— Vv 
1 Y %4 SLT ge N38 oa as 
. 1 
Xp = -—We xX, tyxX%.—pk X3 Xs + Vas 
ep Mie 30 
Be in ek wp 8 (30) 
. =0 = Ms 
x4 = = X2 — (Ws — P X5)X3 
; 1 1 f 
Xs = Flom — Grp 7s 


The expression of the electromagnetic torque becomes as (31). In this work, the goal is to control the speed of 
the motor, in the sense of controlling the flow of the pump. 


3 DMsr 
Rem = =" (xgx3) = KelgsOar = Kemlgs a 
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Based on the reference power and the torque of pump we can deduce the optimum rotation speed of 
the induction motor. The useful power of the induction motor can be written as (32). 


Py = Lemw (32) 
We note 77 the efficiency of the pump motor assembly. It is defined by (33). 


gece (33) 


The torque of pump is expressed as (34) [23]. 
I, = kyw? (34) 


However, the rotational speed can be computed as (35). 


_ 3/9Ppv 
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4.3. Centrifugal pump model 

The centrifugal pumps are widely used for water pumping applications because they are designed for 
high flow-rates and low or medium depths (10 to 100 m). The centrifugal pump is known by the head-flow 
rate H(Q) characteristic curve which depends on the motor speed value as displayed in Figure 11 [24]. The 
flow is proportional to the rotational speed motor while the head manometric (HMT) pressure gauge is 
proportional to the square of the speed. 

The characteristic of the flow head of a centrifugal pump may be approximated by a quadratic function 
using the Pfleider-Peterman model [25], in which the speed of the rotor wis considered as a parameter: 


Hyp = @,W; — ,w,Q — a3Q? (36) 


where a, a2 and a3 are parameters which characterize the pump and given in the technical sheet. The 
characteristic H (Q) of the pipe can be described by (37). 


H=H, + K,Q? (37) 
The HMT (Q) curve as shown in Figure 11, is a parabola. The point where the HMT (Q) curve meets 
the ordinate axis is the point at zero flow. This is the closed valve point or the bubbling point. The operating 


point of the pump is given by the intersection of the characteristic of the water pipe with that of the pump that 
depends on the speed of motor rotation. 
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Figure 11. Characteristic of pump 


5. SIMULATION RESULTS 

The PV Generator used consists of twenty 1Soltech 1STH-215-P type photovoltaic panels. They are 
connected in two parallel strings; each one includes ten panels connected in serial. The boost converter 
parameters are chosen as L = 10 mH, R, = 0.0129,C, = 500 uF and C, = 100 wF.The parameters of a panel 
are given in Table 1. The panel's fuzzy model consists of the following four local models: 
i) G=1000w/m2,T = 5°C,Vypp = 31.71V, Iypp = 7.2A, P = 229.75 W. 
ii) G=1000wW/m2,T = 75°C, Vypp = 22.35V, Iypp = 7.5A, P = 167.8 W. 
ili) G=200w/Am2,T = 5°C,Vypp = 31.62V, Iypp = 1.4554, P = 46 W. 
iv) G=200w/m2,T = 75°C, Vypp = 21.74V, Iypp = 1.513A, P = 32.892 W. 


Table 1. Parameters for 1Soltech 1STH-215-P module 


Parameters at 25 °C and 1000 w/m” for one panel Nominal value 
Maximum power P 213.15 W 
Short circuit current /,.. 784A 
Open circuit voltage V., 36.3 V 
Optimal voltage 29 V 
Optimal current 7.35 A 


The variation of irradiation and temperature are given by the Figures 12 and 13. In the simulation, all 
losses are neglected and it is assumed that the efficiency of the motor pump is equal to unity n=. 
The following Figures 14 to 20 show respectively the evolution of the power delivered by the photovoltaic 
panels, the actual and reference current in the inductance, the flux along the direct axis, the current along the d 
axis, the current along the q axis, the torque and the rotation speed of the induction motor. 

Figures 12 and 13 show that the irradiation varied between 300 to 1000 and the temperature was kept 
constant. It is very logical that the power converted by the panels vary. When irradiation increases the power 
increases and vice versa. Figure 14 clearly appears that the inductance current follows its reference created by the 
supervisor by using the sliding mode control. Figure 16 shows that the field along the direct axis is constant. The 
Figure 17 appears that the current along the direct axis is constant. Figures 18 and 19 show respectively that the 
quadrature current and the torque are proportional. Figure 20 appears that the speed tracks the reference and varies 
smoothly as a function of the power supplied to the motor. Table 2 lists the induction motor parameters. 
Centrifugal pump parameters are: 


w, = 144rad/s,a, = 4.9234 1074 m/(rad/s)? ,az = 1.5826 10-5 m/(rad/s)? (m?/s), 
az = 0.18144 m/(m?/s)’, ky = 0.9,g = 9.81m*/s, ky = 4.4980e-, p = 1000 kg/m? 


Table 2. Induction motor parameters 
Parameters at 25 °C and 1000 w/m? for one panel Nominal value Parameters at 25 °C and 1000 w/ m2 for one panel Nominal value 


Maximum power P 4KW Mutual inductance 0.1722 H 
Stator resistance R, 1.405 Q Moment inertia 0.0131 kg.m? 
Stator inductance L, 5.839 mH Friction 0.002985 rad/s 
Rotor resistance R, 1.395 Q Stator voltage 400 V 

Rotor inductance L, 5.839 mH Frequency 50 Hz 
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6. CONCLUSION 

This paper deals with the optimal operation of a pumping system powered by a photovoltaic source 
using a pump coupled to an induction motor. Conventional and intelligent control techniques have been 
combined to ensure optimal system operation such as sliding mode control, field-oriented control and fuzzy 
systems. Climatic variables such as irradiation and temperature serve as inputs to the fuzzy supervisor which 
in turn provides the optimal values of the photovoltaic source which are current, voltage and power. The MPPT 
algorithm combines T-S fuzzy logic with sliding mode control techniques to extract the maximum energy from 
PV panels. The centrifugal pump is driven by an induction motor, which operates under field-oriented control. 
The flow of the pump varies according to the speed of rotation in a nonlinear way. At each speed of rotation, 
the pump has a well-defined flow rate. We use fuzzy systems to predict the reference speed of the induction 
motor and that of the pump according to the climatic variables temperature and irradiation which in turn make 
it possible to predict the reference power of the PV panels. The simulation results show the effectiveness of 
the proposed approaches. 
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